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Abstract 

We review the calculation for Higgs production via the exclusive re- 
ly-^ ] action pp — ► p + H + p. In the first part we review in some detail 

the calculation of the Durham group and emphasise the main areas of 
; uncertainty. Afterwards, we comment upon other calculations. 

^ ■ 1 Introduction 

Our aim is to compute the cross-section for the process pp — ► p + H + p. We shall only be 
| interested in the kinematic situation where all three final state particles are very far apart in 
rapidity with the Higgs boson the most central. In this "diffractive" situation the scattering 
protons lose only a very small fraction of their energy, but nevertheless enough to produce 
the Higgs boson. Consequently, we are in the limit where the incoming protons have energy 
E much greater than the Higgs mass m# and so we will always neglect terms suppressed by 
powers of mn/E. In the diffractive limit cross-sections do not fall as the beam energy increases 
as a result of gluonic (spin-1) exchanges in the t-channel. 

Given the possibility of instrumenting the LHC to detect protons scattered through tiny 
t£h angles with a high resolution [1-4], diffractive production of any central system X via pp — ► 

p+X+p is immediately of interest if the production rate is large enough. Even if X is as routine 



as a pair of high p T jets we can learn a great deal about QCD in a new regime [2, 3, 5, 6]. But no 
doubt the greatest interest arises if X contains "new physics" [7-19]. The possibility arises to 
measure the new physics in a way that is not possible using the LHC general purpose detectors 
alone. For example, its invariant mass may be measured most accurately, and the spin and CP 
properties of the system may be explored in a manner more akin to methods hitherto thought 



possible only at a future linear collider. Our focus here is on the production of a Standard 
Model Higgs boson [7, 8, 13, 18, 19]. Since the production of the central system X effectively 
factorizes, our calculation will be seen to be of more general utility. 

Most of the time will be spent presenting what we shall call the "Durham Model" of 
central exclusive production [7, 8]. It is based in perturbative QCD and is ultimately to be 
justified a posteriori by checking that there is not a large contribution arising from physics 
below 1 GeV. A little time will also be spent explaining the non-perturbative model presented 
by the Saclay group [13] and inspired by the original paper of Bialas and Landshoff [20]. Even 
less time will be devoted to other approaches which can be viewed, more-or-less, as hybrids of 
the other two [18,19]. 

Apart from the exclusive process we study here, there is also the possibility to produce 
the new physics in conjunction with other centrally produced particles, e.g. pp — > p + H + 
X + p. This more inclusive channel typically has a much higher rate but does not benefit from 
the various advantages of exclusive production. Nevertheless, it must be taken into account 



in any serious phenomenological investigation into the physics potential of central exclusive 
production [21,22] 



2 The Durham Model 

The calculation starts from the easier to compute parton level process qq — > q + H + q shown 
in Figure EJ The Higgs is produced via a top quark loop and a minimum of two gluons need 
to be exchanged in order that no colour be transferred between the incoming and outgoing 
quarks. Quark exchange in the t-channel leads to contributions which are suppressed by an 
inverse power of the beam energy and so the diagram in Figure [T] is the lowest order one. Our 
strategy will be to compute only the imaginary part of the amplitude and we shall make use of 
the Cutkosky rules to do that - the relevant cut is indicated by the vertical dotted line in Figure 
U There is of course a second relevant diagram corresponding to the Higgs being emitted from 
the left-hand gluon. We shall assume that the real part of the amplitude is negligible, as it will be 
in the limit of asymptotically high centre-of-mass energy when the quarks are scattered through 
small angles and the Higgs is produced centrally. 
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Fig. 1: The relevant Feynman graph for qq — > q + H + q. 

The calculation can be further simplified by making use of the eikonal approximation for 
those vertices which couple the gluons to the external quarks. The gluons are very soft and so, 
modulo corrections which are suppressed by the inverse of the beam energy, we can approxi- 
mate the qqg vertices by 2gr^qi t 25x,x>, where r a is a Gell-Mann matrix, g is the QCD coupling 
and the Kronecker delta tells us that the quark does not change its helicity. The calculation of 
the amplitude is now pretty straightforward: 

lmA% = ^x2 Jd(PS) 2 6(( qi -Q) 2 )6((q 2 + Q) 2 ) 
2gq° 2gq 2a 2gqf 2gq^ 

Q2 1.2 1,2 \ui im jn'mk'nl ' V ' 

1 2 

The factor of 1/2 is from the cutting rules and the factor of 2 takes into account that there are 
two diagrams. The phase- space factor is 

d(PS) 2 = S -J dadp (2) 
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where we have introduced the Sudakov variables via Q = aqi + (3q 2 + Qt- The delta functions 
fix the cut quark lines to be on-shell, which means that a —f3 Qt 2 /s <C 1 and Q 2 ~ 
= — Qt 2 - As always, we are neglecting terms which are energy suppressed such as the 
product a(3. For the Higgs production vertex we take the Standard Model result: 

V* - (,„ - V (3, 

where V = m 2 H a s / \Anv)F(m 2 H /m 2 ) and F ss 2/3 provided the Higgs is not too heavy. The 
Durham group also include a NLO K-factor correction to this vertex. After averaging over 
colours we have 

<j ab 

T C T° T a T b — > — 

im jn mk nl a t\j 2 ' 

/ ±iV c 

We can compute the contraction q\Vj£q% either directly or by utilising gauge invariance 
which requires that titVfi = k%V$ = 0. Writing 1 h = xiqi + k iT yields 



v 



n^V ab n u ~ '" 1T ^ 2T \f ab ~ S UP h u V ab <A\ 
Ql V^2 ~ — — V ~ ZT fc lT%V w 

since 2/cx • A; 2 ps a^a^s ~ m 2 H . Note that it is as if the gluons which fuse to produce the Higgs 
are transversely polarized, ~ k iT . Moreover, in the limiting case that the outgoing quarks 
carry no transverse momentum Q T = —k 1T = k 2T and so e x = — e 2 . This is an important 
result; it clearly generalizes to the statement that the centrally produced system should have a 
vanishing ^-component of angular momentum in the limit that the protons scatter through zero 
angle (i.e. qf T <C Q T ). Since we are experimentally interested in very small angle scattering 
this selection rule is effective. One immediate consequence is that the Higgs decay to 6-quarks 
may now be viable. This is because, for massless quarks, the lowest order qq background 
vanishes identically (it does not vanish at NLO). The leading order bb background is therefore 
suppressed by a factor ~ ml/m 2 H . Beyond leading order, one also needs to worry about the bbg 
final state. 

Returning to the task in hand, we can write the colour averaged amplitude as 

^ ~ ^ " 1 x 4a 2 s ( r^;^ V. (5) 



s N 2 S J Q T 2 k 1T 2 k 2T 2 m 2 H 

Using rf 3 qi / rf 3 q 2 / rf 3 q H 5 (4) (gi + q 2 - q[ - q 2 ~ Qh) = d 2 q 1T 'd 2 q 2 T'dy E H (y is the rapidity of 
the Higgs) the cross-section is therefore 



da 




d 2 c LlT 'd 2 c l2T 'dy \ N 2 J (27r)V2 



<^ 2 Qt ' k/JT 2 

2tt Qi ,2 kix 2 k2T 2 3 



(6) 



and for simplicity here we have taken the large top mass limit of V (i.e. m t 3> m H )- We are 
mainly interested in the forward scattering limit whence 



•■IT ' **2T 



1 



Q T 2 ki T 2 k 2 T 2 Qt 4 



It 

'We can do this because Xi ~ mn/y/s whilst the other Sudakov components are ~ Q^/s. 
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As it stands, the integral over Q T diverges. Let us not worry about that for now and instead turn 
our attention to how to convert this parton level cross-section into the hadron level cross-section 
we need. 2 




Fig. 2: The recipe for replacing the quark line (left) by a proton line (right). 

What we really want is the hadronic matrix element which represents the coupling of 
two gluons into a proton, and this is really an off-diagonal parton distribution function [23]. At 
present we don't have much knowledge of these distributions, however we do know the diagonal 
gluon distribution function. Figure |2] illustrates the Durham prescription for coupling the two 
gluons into a proton rather than a quark. The factor K would equal unity if x' = x and kx = 
which is the diagonal limit. That we should, in the amplitude, replace a factor of a s CV/7r by 
dG(x, Qt) /din Q\ can be easily derived starting from the DGLAP equation for evolution off 
an initial quark distribution given by q(x) — 5(1 — x). The Durham approach makes use of a 
result derived in [24] which states that in the case x' <^ x and k\ -C Q% the off-diagonality can 
be approximated by a multiplicative factor, K. Assuming a Gaussian form factor suppression 
for the /c T -dependence they estimate that 

y/ir r(A + 4) U) 

and this result is obtained assuming a simple power-law behaviour of the gluon density, i.e. 
G(x,Q) ~ x~ x . For the production of a 120 GeV Higgs boson at the LHC, K ~ 1.2 x 
e -fefc T /2 j n ^ cross - sec tion, the off-diagonality therefore provides an enhancement of (1.2) 4 ps 
2. Clearly the current lack of knowledge of the off-diagonal gluon is one source of uncertainty 
in the calculation. We also do not really know what to take for the slope parameter b. It should 
perhaps have some dependence upon Qt and for Qt ~ 1.5 GeV, which it will turn out is typical 
for a 120 GeV scalar Higgs, one might anticipate the same /^-dependence as for diffractive J/ ip 
production which is well measured, i.e. 6^4 GeV~ 2 . 

Thus, after integrating over the transverse momenta of the scattered protons we have 



da 1 a s G Fy /2 r '"^ 



dy 256nb 2 9 



f(xi,QT)f(x 2 ,QT) 



(8) 



where f(x, Q) = dG(x, Q)/d\nQ 2 and we have neglected the exchanged transverse momen- 
tum in the integrand. Notice that in determining the total rate we have introduced uncertainty 

2 We note that was first derived by Bialas and Landshoff, except that they made a factor of 2 error in the 
Higgs width to gluons. 
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in the normalisation arising from our lack of knowledge of b. This uncertainty, as we shall soon 
see, is somewhat diminished as the result of a similar 6-dependence in the gap survival factor. 

Now it is time to worry about the fact that our integral diverges in the infra-red. Fortu- 
nately we have missed some crucial physics. The lowest order diagram is not enough, virtual 
graphs possess logarithms in the ratio Qt/^h which are very important as Qt — > 0; these 
logarithms need to be summed to all orders. This is Sudakov physics: thinking in terms of real 
emissions we must be sure to forbid real emissions into the final state. Let's worry about real 
gluon emission off the two gluons which fuse to make the Higgs. The emission probability for 
a single gluon is (assuming for the moment a fixed coupling a s ) 

C A a s r m l/* d$ f m » /2 dE ^ Ca^ ^ /m|\ 
7T JqI v\ J>pt E ~ 4tt 11 \ Q\ ) ' 

The integration limits are kinematic except for the lower limit on the px integral. The fact 
that emissions below Qt are forbidden arises because the gluon not involved in producing the 
Higgs completely screens the colour charge of the fusing gluons if the wavelength of the emitted 
radiation is long enough, i.e. if pt < Qt- Now we see how this helps us solve our infra-red 
problem: as Q T — > so the screening gluon fails to screen and real emission off the fusing 
gluons cannot be suppressed. To see this argument through to its conclusion we realise that 
multiple real emissions exponentiate and so we can write the non-emission probability as 

As Qt — > the exponent diverges and the non-emission probability vanishes faster than any 
power of Q T . In this way our integral over Q T becomes 

J^ff(xi,QT)f(x 2 ,Q T )e- s (10) 

which is finite. 

There are two loose ends to sort out before moving on. Firstly, note that emission off the 
screening gluon is less important since there are no associated logarithms in m# / Qt- Secondly, 
© is correct only so far as the leading double logarithms. It is of considerable practical impor- 
tance to correctly include also the single logarithms. To do this we must re-instate the running 
of a s and allow for the possibility that quarks can be emitted. Including this physics means we 
ought to use 

s ( f m » /4 dPT <*s{Pt) f'- A , r p f t \\\ nn 

e =exp(-/ Q| ^T—l dz[ Z P gg ( z ) + Y:P q9 ( z )}) (11) 

where A = 2p T /m#, and P gg (z) and P qg (z) are the leading order DGLAP splitting functions. 
To correctly sum all single logarithms requires some care in that what we want is the distribution 
of gluons in Q T with no emission up to m H , and this is in fact [25] 

/>, Qt) = -j^yr {e~ S/2 G(x, Q T )) - 
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The integral over Q T is therefore 

J^ff(x u Q T )f(x 2 ,Q T ) (12) 

which reduces to (TTOb in the double logarithmic approximation where the differentiation of the 
Sudakov factor is subleading. 
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Fig. 3: The Higgs cross-section at zero rapidity, and the result obtained if one were to assume that dG(x, Q) /dQ = 
or that dS/dQ = 0. 

The numerical effect of correctly including the single logarithms is large. For production 
of a 120 GeV Higgs at the LHC, there is a factor ~ 30 enhancement compared to the double 
logarithmic approximation, with a large part of this coming from terms involving the deriva- 
tive of the Sudakov. Figure shows just how important it is to keep those single logarithmic 
terms coming from differentiation of the Sudakov factor. For the numerical results we used the 
MRST2001 leading order gluon [26], as included in LHAPDF [27]. Here and elsewhere (unless 
otherwise stated), we use a NLO QCD K-factor of 1.5 and the one-loop running coupling with 
Uf — 4 and Aqcd = 160 MeV. As discussed in the next paragraph, we also formally need an 
infra-red cut-off Q for the QT-integra\; we take Q = 0.3 GeV although as we shall see re- 
sults are insensitive to Q provided it is small enough. Finally, all our results include an overall 
multiplicative "gap survival factor" of 3% (gap survival is discussed shortly). 

Formally there is the problem of the pole in the QCD coupling at p T = Aqcd- However, 
this problem can be side-stepped if the screening gluon has "done its job" sufficiently well and 
rendered an integrand which is peaked at Q T ^> Aqcd since an infra-red cutoff on px can then 
safely be introduced. We must be careful to check whether or not this is the case in processes 
of interest. Indeed, a saddle point estimate of (flOl) reveals that 

exp((lnQ r ))~^exp(--) (13) 
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Fig. 4: The integrand of the Qt integral for three different treatments of a s and mjj = 120 GeV. 

where c is a constant if the gluon density goes like a power of Q T . Clearly there is a tension 
between the Higgs mass, which encourages a large value of the loop momentum, and the sin- 
gular behaviour of the 1 / Q T factor which encourages a low value. Also, as a s reduces so real 
emission is less likely and the Sudakov suppression is less effective in steering Q T away from 
the infra-red. Putting in the numbers one estimates that exp((mQ§,}) ~ 4 GeV 2 for the pro- 
duction of a 120 GeV scalar at the LHC which is just about large enough to permit an analysis 
using perturbative QCD. Figure |4] provides the quantitative support for these statements in the 
case of a Higgs of mass 120 GeV. The integrand of the Qt integral in equation (fT2l is shown 
for both running and fixed a s . We see that the integrand peaks just above 1 GeV and that the 
Sudakov factor becomes increasingly effective in suppressing the cross-section as a s increases. 
Although it isn't too easy to see on this plot, the peak does move to higher values of Qt as 
a s increases in accord with (fT3t . This plot also illustrates quite nicely that the cross-section 
is pretty much insensitive to the infra-red cutoff for Q < 1 GeV and this is made explicit in 
Figure |3 

Discussion of the infra-red sensitivity would not be complete without returning to the 
issue of the unintegrated gluon density. In all our calculations we model the off-diagonality 
as discussed below equation © and we shan't discuss this source of uncertainty any further 
here. 3 Figure |6] shows the gluon density G(x, Q) as determined in four recent global fits (rather 
arbitrarily chosen to illustrate the typical variety) [26, 28-30]. Apart from the Fermi2002 fit, 
they are all leading order fits. Now, none of these parameterisations go down below Q = 1 
GeV, so what is shown in the figure are the gluons extrapolated down to Q — 0. We have 
extrapolated down assuming that the gluon and its derivative are continuous at Q = 1 GeV and 

3 We actually assume a constant enhancement factor of 1.2 per gluon density. 
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Fig. 5: The Higgs cross-section dependence upon the infra-red cutoff Qq. 



that G(x, Q) ~ Q 2 at Q — > 0. 4 The gluons plotted in Figure|6]are all determined at a; = 0.01 
which would be the value probed in the production of a 120 GeV Higgs at y = at the LHC. 
The key point is to note that it is hard to think of any reasonable parameterisation of the gluon 
below 1 GeV which could give a substantial contribution to the cross-section. The Sudakov 
factor suppresses the low Q 2 region and also the size of the gluon and its derivative are crucial, 
and one cannot keep both of these large for Q < 1 GeV. Figure |7] shows the integrand of the 
Qt integral for different fits to the gluon. In all cases the contribution below 1 GeV is small, 
although there are clearly important uncertainties in the cross-section. These uncertainties are 
better seen in Figure [8] which illustrates that one might anticipate a factor of a few uncertainty 
from this source. 

We note that although a variety of parameterizations are presented in Figure [8] the way 
that the actual Qt dependence of the integrand is obtained is the same in each case. In [31, 32] 
the uncertainties arising from the way the unintegrated parton densities are obtained from the 
integrated ones are examined. Here we have followed the prescription presented in [33] which 
amounts to performing one backward step in a DGLAP parton shower. However, it is known 
that such showers tend to underestimate the hardness of, for example, the W/ Z p± spectra in 
hadron colliders unless a large intrinsic transverse momentum is added to the perturbative k± 
distribution of the colliding partons [34]. In [32] it was shown that adding such an intrinsic 
transverse momentum would harden the Qt distribution of the integrand in (fT2l for small Qt 
which in turn lowers the cross-section for central exclusive Higgs production by a factor 2 (for 
a Gaussian intrinsic transverse momentum with (kj_) = 2 GeV 2 ). Investigations into how one 
could use unintegrated gluon densities obtained by CCFM [35] and LDC [36] evolution for 
central exclusive Higgs production have also been performed [32]. However, as discussed in 

4 To be precise we extrapolate assuming G(x, Q) ~ Q 2 +("i r_2 )Q. 
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Fig. 6: The gluon density function in four different parameterisations. 

more detail elsewhere in these proceedings [23], the available parameterizations, which are all 
fitted to HERA data only, are not constrained enough to allow for reliable predictions for Higgs 
production at the LHC. 

This is perhaps a good place to mention pseudo-scalar production, as might occur in an 
extension to the Standard Model. The scalar product, k 1T • k 2 T, in © now becomes (k 1T x 
k 2 T) ■ n, where n is a unit vector along the beam axis. After performing the angular integral the 
only surviving terms are proportional to the vector product of the outgoing proton transverse 
momenta, i.e. qi' x q 2 '. Notice that this term vanishes, in accord with the spin-0 selection rule, 
as q ; ' — > 0. Notice also that the integrand now goes like ~ 1/Qt ( m contrast to the \jQ% in 
the scalar case). As a result c in (fT3l is larger (in fact it is linearly proportional to the power of 
Q T ) and the mean value of Q T smaller. This typically means that pseudo-scalar production is 
not really accessible to a perturbative analysis. 

The Sudakov factor has allowed us to ensure that the exclusive nature of the final state is 
not spoilt by perturbative emission off the hard process. What about non-perturbative particle 
production? The protons can in principle interact quite apart from the perturbative process dis- 
cussed hitherto and this interaction could well lead to the production of additional particles. We 
need to account for the probability that such emission does not occur. Provided the hard process 
leading to the production of the Higgs occurs on a short enough timescale, we might suppose 
that the physics which generates extra particle production factorizes and that its effect can be 
accounted for via an overall factor multiplying the cross-section we have just calculated. This 
is the "gap survival factor". Gap survival is discussed in detail elsewhere in these proceedings 
and so we'll not dwell on it here [37]. 

The gap survival, S 2 , is given by 

da(p + H + p\no soft emission) = da{p + H + p) x S 2 
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Fig. 7: The integrand of the Qt integral for four recent global fits to the gluon. 

where da(p + H + p) is the differential cross-section computed above. The task is to estimate 
S 2 . Clearly this is not straightforward since we cannot utilize QCD perturbation theory. Let us 
at this stage remark that data on a variety of processes observed at HERA, the Tevatron and the 
LHC can help us improve our understanding of "gap survival". 

The model presented here provides a good starting point for understanding the more so- 
phisticated treatments [38-40]. Dynamically, one expects that the likelihood of extra particle 
production will be greater if the incoming protons collide at small transverse separation com- 
pared to collisions at larger separations. The simplest model which is capable of capturing this 
feature is one which additionally assumes that there is a single soft particle production mech- 
anism, let us call it a "re- scattering event", and that re-scattering events are independent of 
each other for a collision between two protons at transverse separation r. In such a model we 
can use Poisson statistics to model the distribution in the number of re-scattering events per 
proton-proton interaction: 

P„(r) = ^exp(- X (r)). (14) 
n! 

This is the probability of having n re-scattering events where %(r) is the mean number of such 
events for proton-proton collisions at transverse separation r. Clearly the important dynamics 
resides in x( r ); we expect it to fall monatonically as r increases and that it should be much 
smaller than unity for r much greater than the QCD radius of the proton. Let us for the moment 
assume we know x( r )> then we can determine S* 2 via 

S 2 = I dr d(T ( r ) expt-xO)) (15) 
/ dr da(r) 

where dcr(r) is the cross-section for the hard process that produces the Higgs expressed in 
terms of the transverse separation of the protons. Everything except the r dependence of da 
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Fig. 8: The Higgs cross-section for four recent global fits to the gluon. 

cancels when computing S 2 and so we need focus only on the dependence of the hard process 
on the transverse momenta of the scattered protons (qi'), these being Fourier conjugate to the 
transverse position of the protons, i.e. 

da(r) oc [(J riV e iqiW2 exp(-6 qi /2 /2)) x ( / ri 2 q 2 ' e^'^ 2 exp(-6q 2 /2 /2))] 2 

oc ex P (- r -^j . (16) 

Notice that since the b here is the same as that which enters into the denominator of the expres- 
sion for the total rate there is the aforementioned reduced sensitivity to b since as b decreases so 
does S 2 (since the collisions are necessarily more central) and what matters is the ratio S 2 /b 2 . 

It remains for us to determine the mean multiplicity x(r). If there really is only one type 
of re-scattering event 5 independent of the hard scattering, then the inelastic scattering cross- 
section can be written 

^inelastic = / ^M 1 ~ eX P (~ X( r ))) , O 7 ) 

from which it follows that the elastic and total cross-sections are 

Elastic = / d 2 r(l - exp(- X (r)/2)) 2 , (18) 
fftotai = 2y>r(l-exp(- X (r)/2)). (19) 

There is an abundance of data which we can use to test this model and we can proceed to perform 
a parametric fit to x( r )- This is essentially what is done in the literature, sometimes going 

5 Clearly this is not actually the case, but such a "single channel eikonal" model has the benefit of being simple. 
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beyond a single-channel approach. Suffice to say that this simple approach works rather well. 
Moreover, it also underpins the models of the underlying event currently implemented in the 
PYTHIA [41] and HERWIG [42,43] Monte Carlo event generators which have so far been quite 
successful in describing many of the features of the underlying event [44-46]. Typically, models 
of gap survival predict S 2 of a few percent at the LHC. Although data support the existing 
models of gap survival there is considerable room for improvement in testing them further and 
in so doing gaining greater control of what is perhaps the major theoretical uncertainty in the 
computation of exclusive Higgs production. In all our plots we took S 2 = 3% which is typical 
of the estimates in the literature for Higgs production ath the LHC. 

3 Other Models 

We'll focus in this section mainly on the model presented by what we shall call the Saclay 
group [13]. The model is a direct implementation of the original Bialas-Landshoff (BL) cal- 
culation [20] supplemented with a gap survival factor. It must be emphasised that BL did not 
claim to have computed for an exclusive process, indeed they were careful to state that "ad- 
ditional. ..interactions. ..will generate extra particles. ..Thus our calculation really is an inclusive 
one". 

Equation © is the last equation that is common to both models. BL account for the 
coupling to the proton in a very simple manner: they multiply the quark level amplitude by 
a factor of 9 (which corresponds to assuming that there are three quarks in each proton that 
are able to scatter off each other). Exactly like the Durham group they also include a form 
factor suppression factor exp(— bq'^) for each proton at the cross-section level with 6 = 4 
GeV 2 . Since BL are not interested in suppressing radiation, they do have a problem with 
the infra-red since there is no Sudakov factor. They dealt with this by following the earlier 
efforts of Landshoff and Nachtmann (LN) in replacing the perturbative gluon propagators with 
non-perturbative ones [47,48]: 

Rather arbitrarily, g 2 = An was assumed, except for the coupling of the gluons to the top quark 
loop, where a s — 0.1 was used. 

Following LN, /i and A are determined by assuming that the pp elastic scattering cross- 
section at high energy can be approximated by the exchange of two of these non-perturbative 
gluons between the 3x3 constituent quarks: the imaginary part of this amplitude determines 
the total cross-section for which there are data which can be fitted to. In order to carry out this 
procedure successfully, one needs to recognize that a two-gluon exchange model is never going 
to yield the gentle rise with increasing centre-of-mass energy characteristic of the total cross- 
section. BL therefore also include an additional "reggeization" factor of s^*) -1 in the elastic 
scattering amplitude where 

a(t) = 1 + e + a' t 

is the pomeron trajectory which ensures that a good fit to total cross-section data is possible 
for e = 0.08 and a' = 0.25 GeV~ 2 . In this way the two-gluon system is modelling pomeron 
exchange. They found that [i w 1 GeV and A « 30 GeV 2 gave a good fit to the data. Similarly, 
the amplitude for central Higgs production picks up two reggeization factors. 
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The inclusive production of a Higgs boson in association with two final state protons is 
clearly much more infra-red sensitive than the exclusive case where the Sudakov factor saves 
the day. Having said that, the Saclay model does not include the Sudakov suppression factor. 
Instead it relies upon the behaviour of the non-perturative gluon propagators to render the Qt 
integral finite. As a result, the typical Qt is much smaller than in the Durham case. Indeed it 
may be sufficiently small to make the approximation Q\ ^> q^ T invalid which means that the 
spin-0 selection rule is no longer applicable. 

Pulling everything together, the Saclay model of the cross-section for pp — > p + H + p 

gives 



da 



d 2 q 1T 'd 2 q2 T 'dy 



S 2 



iv c 2 -i 

N 2 
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7T 



*1 x 2 



d 2 QT n 2 (AY _ 2/ 2 ,2 



(20) 



The reggeization factors depend upon the momentum fractions x\ and x 2 which satisfy x\X 2 s = 
m 2 H and y = \ hi(x\/x2). The only difference 6 between this and the original BL result is the 
factor of S* 2 . Integrating over the final state transverse momenta and simplifying a little gives 
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A 3 fi 4 
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Fig. 9: Comparing dependence upon ran of the Saclay and Durham predictions. S 2 = 3% in both cases. 



Apart from the factor 2 error previously mentioned. 
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Fig. 10: Comparing dependence upon y/s of the Saclay and Durham predictions for mjj = 120 GeV. 

Figure |9] shows how the Saclay model typically predicts a rather larger cross-section with 
a weaker dependence upon m# than the Durham model. The weaker dependence upon m H 
arises because the Saclay model does not contain the Sudakov suppression, which is more 
pronounced at larger tuh, and also because of the choice e = 0.08. A larger value would induce 
a correspondingly more rapid fall. The Durham use of the gluon density function does indeed 
translate into an effective value of e subtantially larger than 0.08. This effect is also to be seen 
in the dependence of the model predictions upon the centre-of-mass energy as shown in Figure 
ITOl We have once again assumed a constant S 2 = 3% in this figure despite the fact that one 
does expect a dependence of the gap survival factor upon the energy. 

Figure [TT] compares the rapidity dependence of the Higgs production cross-section in the 
two models. The Saclay prediction is almost y-independent. Indeed the only y-dependence is 
a consequence of a' ^ 0. In both models the calculations are really only meant to be used for 
centrally produced Higgs bosons, i.e. \y \ not too large since otherwise one ought to revisit the 
approximations implicit in taking the high-energy limit. Nevertheless, the Durham prediction 
does anticipate a fall as \y\ increases, and this is coming because one is probing larger values 
of x in the gluon density. In contrast, the Saclay prediction does not anticipate this fall and so a 
cutoff in rapidity needs to be introduced in quoting any cross-section integrated over rapidity. In 
FigurefTTIa cut on 2i 2 < 0.1 is made (which is equivalent to a cut on \y\ < 2.5) for the Saclay 
model. After integrating over rapidity, the Durham model predicts a total cross-section of 2 fb 
for the production of a 120 GeV Higgs boson at the LHC whilst the Saclay model anticipates a 
cross-section a factor ~ 5 larger. 

The essentially non-perturbative Saclay prediction clearly has some very substantial un- 
certainties associated with it. The choice of an exponentially falling gluon propagator means 
that there is no place for a perturbative component. However, as the Durham calculation shows, 
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Fig. 11: Comparing the y dependence of the Saclay and Durham predictions for run = 120 GeV. Note that the 
Saclay prediction has been reduced by a factor 5 to make the plot easier to read. The numbers in parenthesis are 
the total cross-sections, i.e. integrated over rapidity. 

there does not seem to be any good reason for neglecting contributions from perturbatively 
large values of Q T . It also seems entirely reasonable to object on the grounds that one should 
not neglect the Sudakov suppression factor and that including it would substantially reduce the 
cross-section. 

In [18], the Sudakov factor of equation dTTb is included, with the rest of the amplitude 
computed following Bialas-Landshoff. The perturbative Sudakov factor is also included in the 
approach of [19], albeit only at the level of the double logarithms. This latter approach uses 
perturbative gluons throughout the calculation but Regge factors are included to determine the 
coupling of the gluons into the protons, i.e. rather than the unintegrated partons of the Durham 
model. In both cases the perturbative Sudakov factor, not suprisingly, is important. 

4 Concluding remarks 

We hope to have provided a detailed introduction to the Durham model for central exclusive 
Higgs production. The underlying theory has been explained and the various sources of un- 
certainty highlighted with particular emphasis on the sensitivity of the predictions to gluon 
dynamics in the infra-red region. We also made some attempt to mention other approaches 
which can be found in the literature. 

The focus has been on the production of a Standard Model Higgs boson but it should 
be clear that the formalism can readily be applied to the central production of any system X 
which has a coupling to gluons and invariant mass much smaller than the beam energy. There 
are many very interesting possibilities for system X which have been explored in the literature 
and we have not made any attempt to explore them here [2, 3, 8, 11, 15-17]. Nor have we 
paid any attention to the crucial challenge of separating signal events from background [5,9]. 
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The inclusion of theoretical models into Monte Carlo event generators and a discussion of the 
experimental issues relating to central exclusive particle production have not been considered 
here but can be found in other contributions to these proceedings [49,50]. 

It seems that perturbative QCD can be used to compute cross-sections for processes of the 
type pp — > p + X + p. The calculations are uncertain but indicate that rates ought to be high 
enough to be interesting at the LHC. In the case that the system X is a pair of jets there ought 
to be the possibility to explore this physics at the Tevatron [51]. Information gained from such 
an analysis would help pin down theoretical uncertainties, as would information on the rarer 
but cleaner channel where X is a pair of photons [52]. Of greatest interest is when X contains 
"new physics" whence this central exclusive production mechanism offers new possibilities for 
its exploration. 
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